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Abstract 

We explore local consequences of a non-zero cosmological constant on astrophysical structures. 
We find that the effects are not only sensitive to the density of the configurations but also to the 
. geometry. For non- homogeneous configurations, we calculate the effects for a polytropic configura- 

tions and the isothermal sphere. Special emphasis is put on the fact that the cosmological constant 
£ \ sets certain scales of length, time, mass and density. Sizable effects are established for non spherical 

systems such as elliptical galaxy clusters where the effects of A are growing with the flatness of the 
s ! ' system. The equilibrium of rotating ellipsoids is modified and the cosmological constant allows new 

^ ' configurations of equilibrium. 

(N 

1 Introduction: cosmological motivations for Dark Energy 

' Astronomical data based on light curves of distant la Supernova £P, anisotropics in the cosmic back- 
, ground radiation [5] and the matter power spectrum of large scale structure [3] agrees with an accelerated 
universe which is spatially fiat and dominated by a Dark Energy component with a contribution of 70%. 
The remaining 30% corresponds to a cold dark matter component, while the contribution of radiation 
and baryons is negligible. In order to account for an accelerated phase at the present epoch, the dark 
energy component must have a negative equation of state p x = uj x p x , uj x < 0. The most favored scenario 
is the model with a cosmological constant A giving p x = p vac = —p V3 , c = A/8ir (we set the Newtonian 
constant, Gn, to one), called ACDM model. However there are other interesting models introduced 
in order to mimic the current contribution of the Dark Energy component. Some of the most relevant 
models are i) Dark Energy with p x — ui K p x where — 1 < lu x < 0, ii) Dynamical Dark Energy (DDE) 
where the equation of state is a function of time ui x = tu x {t), iii) Chaplygin Gas 0] with an equation of 



state of the form p = —czHq ' p^ 1 where a < (3/87r) 7+1 , a consequence of a > and finite age of 
the universe, iv) higher order gravity [5] and v) Quintessence models jS] which introduce the concept 
of scalar fields ruled by a potential V(4>) which at the present time mimics a cosmological constant as 
A V(<po) where tpo is the value of the field at the minimum of the potential. Although in principle 
one could translate one model into the other, it is worth noting that Quintessence models have been not 
only successful when trying to interpret the cosmological constant but also in explaining the well known 
problems of standard cosmology (i.e, the horizon problem, the homogeneity problem and fine-tunning 
problems) [7]. 

After the discovery of the current state of the universe ( the fact that we are dominated by the cosmo- 
logical constant), several question have arisen. The inferred values for cosmological parameters such as 
the density parameter associated with the cold dark matter r2 c a m = 0.27 ± 0.04 and the dark energy 
component associated with e.g. the cosmological constant i7 V ac = 0.73 ± 0.04 lead to the coincidence 
problem. One such coincidence emerges if we compare the scales of length and time for the universe 
with values set by the cosmological constant. The length (time) scale set by A is defined as r\ = A -1 / 2 
which corresponds to 4 x 10 3 Mpc (9.6Gyr). These scales coincide with the size (age) of the universe 
which is proportional to the Hubble's distance (time) cIh = Hq 1 w 6 x 10 3 Mpc (14Gyr). Essentially the 
coincidence is that we are living in an epoch dominated by A where e.g. r\ ~ Hq 1 , a relation possible 
only at certain times of the expansion. Other aspects of some astrophysical scales set by A can be found 
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in [S] . Nowadays the origin of the cosmological constant and related problems are a very active field of 
modern physics UJ. 

In standard cosmology the universe is considered an isotropic and homogeneous fluid at large scales. This 
allows us to describe it by the Robertson-Friedmann- Walker metric ds 2 = —dt 2 +a(t) 2 ((1 — kr 2 )~ 1 dr 2 + r 2 dQ 2 ) 
where k represents the curvature, dfl 2 is the metric on a 2-sphere and a(t) is the dimcnsionlcss scale 
factor. This factor determines the kinematics of the universe though the Einstein field equations. For 
a spatially flat universe one obtains the Friedmann equation (a/a) 2 = H 2 (z) = H^VL vac hi(z) and accel- 
eration equation a/a = HQfl vac h2(z). The time variable z is the red-shift defined as z — a -1 — 1. The 
functions h\ t 2(z) are given as 



hx{z) 



^cdn 



■{l + zf + (l + z) f ^, 



h 2 (z) 



n 



cdm 



(l + z) 3 + (l + z) /(z) (l + 3w x (z)) 



, (1) 



We have only considered contributions from a cold dark matter component and vacuum energy density. 
The function f(z) is given by 



l+w x (z') 



ln(l + z) J 1 + z 



(2) 



—3H(t)(p + p) for the dark energy compo- 
-1 corresponds to the cosmological constant 
,. Alternative models with u> K < — 1 lead to future singularities in the so-called Phantom regime 



which follows from the energy-density conservation law p ■■ 
nent p K (z) = u> K (z)p x (z). As pointed before, the case tu x = 
P*. = 



2 Local Dynamics with Dark Energy 

Local effects of Dark Energy have been investigated by many authors ^J. We explore these effects 
of the background through the Newtonian limit of Einstein field equations |14j from which one can 
derive a modified Poisson's equation V 2 $ = 4ir5p — 3Hoil vac h2(z), where 5p is the ovcrdcnsity that 
gives rise to the gravitational source of the potential The total energy density within the clustered 
configuration is a contribution of the cold dark matter in the background and the collapsed fraction, 
i.e, p(z) = /O c dm(0)(l + z) 3 + Sp. The extra-term 3HaQ vac h2(z) in Poisson's equation reduces to 87T/j V ac 
for uj x = — 1 if we neglect the contribution from the cold dark matter component. This is the standard 
Newtonian limit. The more general limit retaining also the time (z) dependent parts is derived e.g. in 

na. 

The solution for the potential can be simplified as $(r, z) = $gr av (r) — H'^VL vac h2{z)r 2 which defines the 
Newton- Hooke spacetime |13| for u> x = — 1 and p c dm ~ at the present time. Two effects combine in 
order that we are allowed to make this approximation at the present time: i)the red-shift-dependence 
(1 + z) 3 and the ratio p c dm/ P with p the mean density of the configuration. This ratio was 10 -2 at the 
moment of virialisation. 

Once one knows the solution of Poisson equation $(r, z), one can write the Euler's equation for the 
ovcrdcnsity as Spiii + djVij + SpdiQ = 0, (neglecting p c dm), where Vij is the dispersion tensor. Taking 
spatial moments on Euler's equation one can derive the second order tensor virial equation |15| as 



1 d 2 l ik 1 2 . 

- -H a Q vac (l + z) 



2 d* 2 



h 1 (z)(l + z)^ + (h2(z) + 2h 1 (z)) d 



dz 2 v ^ w iv " dz 
= 2T lk -\Wfr\+ H o^Mz)2 lk +Tl lkl (3) 



where we have neglected surface integrals (this is not the case when we consider configurations with a 
non zero pressure at the boundary or if we consider non radial external forces, as electromagnetic forces). 
The first line of is the second derivative with respect to time of the inertial tensor in dependence of 
the red-shift where we have used the Friedmann equation (h\(z)) and the acceleration equation (/i2(z)). 
The second line of iJSJ is the dynamical information encoded in the virial equation. The quantities 
involved in this expression are defined as 1 

l lk = [ pnr k d 3 r, T lk = \ [ pu t u k d 3 r Wf fe rav = - / pnd^^r, U ik = [ T lk d 3 r, (4) 
J_v j Jv Jv Jv 

^^We have renamed the density of the configurations from 5p to p. 
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corresponds to the moment of inertia tensor, the kinetic energy tensor, the gravitational potential energy 
tensor and the dispersion tensor. 

By taking the trace on Eq. J2Jl and assuming the virial graviational equilibrium we obtain the most 
known form for the virial theorem: 

|>V srav | = 2T + 3n + |7rp vac X, (5) 

where T = Tr(Tjfc) and so on for the other terms. Note then that the effects of a positive cosmological 
constant are coupled to the moment of inertia which depends on the density profile (or the equation 
of state) and the geometry of the configuration. Hence, these two parameters (density and geometry) 
become relevant in order to determine sizable effects of a vacuum energy density. It is to be expected 
that some effects due to A will show up only if the density of the object is low. However, we will see 
that geometry helps to enhance such effects. 

In order to describe a general consequence of Eq.JSJ, we write the total kinetic energy as K, = T + 3H. 
Since this quantity is positive defined, we conclude that the state of equilibrium is characterized by the 
inequality g > Ap V3jC where g is a characteristic parameter with units of density (such as the man density 
or the central density) and A = (167r/3)f/|>V grav |, together with 2 = 1/g and |WS rav | = 2\W\/g 2 . The 
parameter A depends on the geometry of the configurations and/or the parameters specifying internal 
structure such as the equation of state. 



3 Effects on astrophysical configurations 

3.1 Spherical configurations 

The tensor virial equation is trivially satisfied for spherically symmetric configuration (in the absence 
of electromagnetic fields or non radial forces) since Wf™ v = <5jfcW srav and 2^. = Therefore Eq.JSJ 

is sufficient in order describe configurations with this geometry. 

3.1.1 Homogeneous systems 

For constant density (with g = p) it is straightforward to check that Asph = 2. In this case the effects 
of Pvac arc not enhanced because of this geometrical factor. However, other effects can be shown to 
exist. Indeed, the presence of a cosmological constant allows the existence of a maximal virial radius of 
a spherical configuration in virial equilibrium in the limit when K, — > 0. The radius at virial equilibrium 
can be calculated from Eq. (JSJ), as the solution of the cubic equation i?^ ir + (lO^r^J i? v i r — 3Mr 2 = 
where 77 = JC/M = (3fcs//x)T, \x is the mass of the average member of the configuration, ks is the 
Boltzmann constant, T is the temperature. The solution for the radius at virial equilibrium reads as 

Rvir{x,rf) = 2.53a;~ 1/3 ?7 1/2 sinh 

where x = Mr^ 1 and i? max = i? vir (a;,0) = {3Mr\) 1 /' i is the maximum radius allowed for astrophysi- 
cal/cosmological structures. In astrophysical terms we have 

i?max - 9.5 X 1(T 2 I— fiva!/l 70 3 k P C. (7) 

This length scale is a combination of a cosmological length scale defined by the cosmological constant 
ta and a pure astrophysical scale, i.e, the Schwarszchild's radius R s = 2M. This combination gives rise 
to a relevant astrophysical scale. This new scale with a replacement M — > \x in (JJJ also appears when 
studying the motion of test particles in the Einstein-de Sitter space time |I6| and represents the last 
bound orbit around the source of gravity with mass p. In other words, Eq. {JJ would imply that the 
cosmological constant sets the typical scale for bound configurations. Some examples are welcome: for 
/.t = I0 6 A/q i.e. a typical mass of a globular clusters (the tight connection between galaxies and globular 
clusters is that the latter are important in galaxy formation) we get i? max ~ 8 kpc, whereas for a galaxy 
mass, fi = lO n M0, we obtain i? max ~3x 10 2 kpc. The numerical outcomes correspond to the radius a 
typical galaxy in the first example, while the second one agrees with numerical values for the extension 
of a galaxy cluster. 



f 
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3.1.2 Non-homogeneous systems 

The effects on non-homogcncous spherical configurations have been explored by considering polytropic 
configurations characterized by a barotropic equation of state of the form p = np 1+1 / n The 
equation governing the behavior of the density of polytropic configurations with cosmological constant 
is the modified Lanc-Emden equation written as 

(<■*)-'-*•■ (8) 

where £ = r/a, ip = (p/p c ) 1/n , p c = p(r = 0), a = («(n + l)/47r) 1 / 2 pc / ' 2n ~ 1 ' /2 and Cc = 2p vac /p c . As a 
main effect of a positive A one can demonstrate that not all the configurations with a polytropic equation 
of state have a defined radius, even for n < 5, as can be seen in Fig^i.c. not all values of n leads to a 
solution ijj such that for some £i one gets ^(£1 = 0)). Wc generalize the equilibrium criteria given for 
homogeneous configurations by identifying g = p c and writing p c > AnPv&c where the factor A n takes 
different values for different indexes n. Some examples are A\ = 10.24, A3/2 = 24.24, A3 = 307.69, 
Aa ?« 4000. Characteristic quantities and concepts such as the radius and the mass as well as the 
stability criteria get modified in the presence of A. For instance, the radius of such configurations can 
be written as R n = a n Ro where Rq is the radius that one would obtain with p vac = and a n = a ra (Cc) 
is an enhancement factor which depends on the polytropic index n and the ratio £ c - Some numerical 
values are ai(0.1) = 1.12, a 3/2 (0.05) = 1.11, ai (0.001) = 1.001, a 3 (0.001) = 1.022 which shows that 
even for spherical configurations effects are found for low density configurations with different values of 
n. Worth noting is the case of the isothermal sphere i.e. the case n — » 00. The isothermal sphere is an 
often employed model in astrophysics ^7]. From Figure 1 we can see that including A the cases n > 5 
do not have a defined radius even in the asymptotic sense. Therefore one can easily suspect that the 
extreme case of the isothermal sphere will display even bigger effects which is to say, sizable effects for 
larger density. This is indeed the case and the isothermal sphere is not a viable model (with A) even at 
density comparable with galaxy densities |19j . The stability criteria for polytropic configurations can 
be derived via the virial theorem in two ways: one option is to develop Lagrangian perturbations of 
the form £e luJt (where £ is the Lagrangian displacement) on the tensor virial equation and derive an 
equation for the oscillation frequency around equilibrium uj, such that stability would implies lo 2 > 
(see for instance Another option is to solve for the mass of the configuration as a function of the 

central density. Stability (instability) stands for dM/dp c > (dM/dp c < 0) |21|. Both ways lead us to 
write a critical value for the polytropic index such that polytropic configuration are stable if 7 > 7 cr it 
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where 



7crit 



4 2 ding 

7crit(Cc) = - + r- 



3 3 <9 In p c ' 



(9) 



is the critical polytropic index |19) . The function Q depends on the ratio £ c through integrals of the 
Lane-Emden function [15] which in turn depends on the index n. Hence, Eq © is a transcendental 
equation for the critical value of the polytropic index. Nevertheless, Eq. © is a generalization of stability 
criteria for homogeneous configurations under radial adiabatic perturbations with cosmological constant 
and can be written as 7 > (4/3)(l + (l/2)£). It is worth mentioning that by including the corrections 
due to general relativity, the critical value for 7 cr ;t is also modified as 7 cr it = (4/3) + R s /R [22] and 
hence for compact objects the correction to the critical polytropic index is stronger from the effects of 
general relativity than from the effects of the background, as expected. 

3.2 Non-spherical configurations 

As mentioned before, the effects of a positive cosmological constant can be enhanced through the geo- 
metrical factor A. As an example, we consider the virial theorem and solve for the velocity dispersion 
for arbitrary geometry and constant density. The final expression has the form v 2 /v\_ = 1 — (\/2)QA 
where C = 2p vac /p. In the last section we saw that for spherical and homogeneous systems -4 sp h = 2 and 
hence the correction is just 1 — C which implies that a sizable effect could be reached only in systems with 
p — > 2p vac . On the other hand, for non-spherical configurations we can combine a low density feature 
with a highly non-spherical geometry in order to find sizable effects on realistic systems. For instance, 
for oblate and prolate homogeneous configurations these geometrical factor is written as 117) 
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(10) 



where e 2 = 1 — a§/a 2 is the eccentricity for oblate configurations and e 2 = 1 — a 2 /a 2 , the corresponding 
one for the prolate case. From Ea. (|10|) we see that very flat astrophysical objects we can gain several 
orders of magnitude of enhancement of the effect of p vac - Clearly we often encounter in the universe flat 
configurations such as elliptical galaxies, spiral disk galaxies, clusters of galaxies of different forms and 
finally superclusters which can have the forms of pancakes. Of course, the more dilute the system is, 
the bigger the effect of A. We can expect sizable effects for clusters and superclusters, even for very flat 
galaxies. Regarding the latter, low density galaxies like the nearly invisible galaxies are among other 
the best candidates. 

In order to explore the effects of p vac further we turn to the tensor virial equation for a rotating oblate 
configuration. Writing IT^ = fl 2 ot Zik — Q, m tiIkjQjmt as (twice) the rotational kinetic energy tensor 
and defining a coordinate system such that (l = Vte z one gets the following expression for the angular 
velocity: 



1 win 



1 win 
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(ii) 



where f2 corresponding to the angular velocity when A 
and C = 2p vac /p. The function g(e) has been defined as 



given by the Maclaurin formula (see |18j ') 



5(e) = ^ 



(l-e 2 ) 1 /2(3-2e 2 ) 



arcsm e 



3e(l - e 2 ) 



(12) 



As can be seen from these expressions, the vacuum energy density p V ac has a two effects on the angular 
velocity. In first case, it reduces the angular velocity with respect to the value VLq especially at the local 
maximum (see Figure l3~2f . This is not a small effect and can affect even galaxies. Secondly, we see from 
Eq. \i2\ that in the limit e — > 1 one has (\/2)Qg(e) — * (l/2)C(327r/9) (03/01), approaching 1 for a very 
flat oblate configuration and not too dense matter. Therefore, beyond the local maximum in tt ro t the 
cosmological constant causes a steeper fall of f2 ro t towards 0. In passing we note that galaxy clusters 
exhibit rotating ellipsoidal configurations j23| to which our results can be applied. 
Figure 2 shows the angular velocity and the ratio j3 defined as the ratio of rotational over gravitational 
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contributions i.e. (3 = T/|W grav | versus the eccentricity e. 

The tensor virial equation not only gives us the equation for angular velocity (which can be also derived 
from Euler's equation) but also constrains the non-spherical geometries in equilibrium. That is, not all 
values of eccentricities represent configurations with rotation which are in equilibrium |20| . For instance, 
with p vac = only two configurations are allowed to have angular velocity along the minor axis and 
to be in equilibrium, namely, the Maclaurin solution (oblate) with < (73 < 1 and the Jacobi solution 
(triaxial) with < (73 < 0.582 where 03 = 03/01 (the value gJp = 0.582 is known as the bifurcation 
point). With the inclusion of p vac in the virial equation, we find two effects |24| . In first place, we see 
that the bifurcation point becomes a increasing function of Secondly, we find a second bifurcation 
point 93(n lin ) such that for q 3 < q}jf mla \ , the only allowed solution is again the Maclaurin solution. 
Similar results are found for triaxial configurations which can be reduced to prolate ellipsoids and angular 
velocity in the direction of the largest axis which is called minor axis rotation. This is rare case, but 
can nevertheless be found in nature |25| . 



4 Rotating configurations in an expanding background 



In this section we consider the full contribution of the background. This is to say we examine how 
virializcd matter responses to a time dependent external force representing the expansion. We assume 
that the shape of the system remains constant even under the effects of a time dependent background. It 
is clear that a time dependent term in the virial equation implies that certain internal properties of the 
matter distribution (like angular velocity and internal velocity) also change in time, in such a way as to 
maintain a constant volume and density. We concentrate here on the effects on angular velocity. Using 
equation (1) and (3) and assuming the left hand side of (3) to be zero (i.e. we assume gravitational 
equilibrium) we obtain 



e(z;w x ,C,e) 



n 2 -n 2 

"0 



^cdr 



a 



■(I + zf + (l + 3Lu x {z)){l + z y^ 



(13) 



where fio is the angular velocity calculated without any Dark Energy background. As mentioned above 
in order to maintain a constant shape in a expanding background the system has to pay the price that 
its angular velocity will change with time, according to the dominant component of the background. 

If we neglect the CDM component with respect to the density of the system, we see that the Dark 
Energy component makes the angular velocity to decrease in order for the system to maintain equilibrium 
since 1 + 3w x < 0. At z = we get e ~ —0Ag(e)(. An extreme situation would be reached if e — ► 1 so 
that O 2 -> 0. For e - 0.8 this would require £ ~ 0.5, which is a very diluted configuration. For higher 
eccentricities (say e ~ 0.97) one would need a system with £ ~ 0.2 which is still a rather low value. In 
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Figure 3: The fractional variation of the angular velocity with respect its value for A = as a function of the rcdshift 
for different Dark Energy models with typical values of £ and eccentricities. 



1/3 

the ACDM model SI = is reached at a redshift given by z c = [(2Jl vac (£g(e) — 2))/(£g(e)Jl c dm)] — 1- 
For realistic examples as £ ~ 0.023 (this represents the case of a configurations with p ~ 200p c d m which 
is the typical value predicted by the nonlinear collapse model) an allowed redshift could be reached if 
e is very close to 1. This clearly means the extreme situation of z c requires extreme flat objects. In 
Fig. 21 we show the behavior of e as a function of the redshift with ( = 0.01,0.1 and e = 0.9,0.95 for 
three different values of the equation of state for Dark Energy. It is clear that the effects associated 
with a cosmological constant are stronger than the ones associated with other Dark Energy models, 
but in general those effects are small for realistic values of £ as the used in the figure. We would have 
to measure the angular velocity at different red-shifts very exactly to see an effect over a range of z. 
However, the difference between the models is more significant. 



5 Conclusions 

In the present work we investigated in detail the astrophysical relevance of the cosmological constant for 
equilibrium configurations. Using the tensor and scalar virial equations and the Lane-Emden equation 
we could show that many astrophysical facets get modified by A. An important aspect is the fact that 
A introduces new relevant scales (these scales would be zero if A = 0) like the maximal virial volume 
defined by the maximal virial radius © and the maximal extension of bound orbits. 

It is often assumed that supcrclustcrs with densities ~ p cr ;t are not in equilibrium. With the 
inequality g > *4p vac we have a precise tool to quantify this statement. Indeed, the pancake structure 
of the supcrclustcrs lead us to the conclusion that they are even far away from the equilibrium state 
due to the factor A which grows with the object's flatness. On the other hand we can apply the results 
also to low density clusters which are assumed to be virialized. Depending on shape and density various 
observable, like internal velocity or angular velocity, will get affected by the repulsive force induced 
by Dark Energy. The effects are bigger the large the geometry of the objects deviates from spherical 
symmetry and the lower the density of the object. Finally, the response of the objects to epoch dependent 
expansion gets reflected also in a time dependent variation of some of their properties and the difference 
between different Dark Energy models becomes manifest. 
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